On the electrical double layer contribution to the interfacial tension of protein crystals 
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We study the electrical double layer at the interface between a protein crystal and a salt solution 
or a dilute solution of protein, and estimate the double layer's contribution to the interfacial tension 
of this interface. This contribution is negative and decreases in magnitude with increasing salt 
concentration. We also consider briefly the interaction between a pair of protein surfaces. 



I. INTRODUCTION 

Protein crystallisation is not only of immense practical 
importance, but also attracts interest as a phenomenon 
in its own right. The conditions under which crystalli- 
sation readily occurs are elusive and vary unpredictably, 
and a number of hypotheses to explain this have been put 
forward |Q, g, Q [| . The process of crystallisation, as 
with other first-order phase transitions, begins with nu- 
cleation, where a microscopic nucleus of the crystalline 
phase first forms. The elementary theory of nucleation, 
generally called classical nucleation theory [Q, predicts 
a rate of nucleation which varies as the exponential of 
minus the cube of the interfacial tension and hence is 
very sensitive to the magnitude of the interfacial tension. 
In the present paper, we extend earlier work of one of 
us on the effect of salt on the bulk phase behaviour 
to calculate the effect of salt on the interfacial tension 
and electrical structure at the surface of a protein crys- 
tal. Our calculation can be compared to the recent work 
by Haas and Drenth who also develop a theory for the 
interfacial tension of protein crystals || . Their theory is 
of the Cahn-Hilliard type and does not treat electrostatic 
effects explicitly 

It is a general consequence of the long range nature of 
the Coulomb law that all bulk phases have to be elec- 
trically neutral and space charge effects are confined to 
interfacial regions. The consequences of charge neutral- 
ity for membrane equilibrium were explored by Frederick 
Donnan as long ago as 1911 1 1 . Much more recently, 



the electrical structure of the interface in an electrolyte 
near the critical point was studied by Nabutovskii and 
Nemov |ll[. In a crystal of charged protein molecules, 
the counterion concentration in the interior of the crystal 
will be higher than in the salt solution outside the crys- 
tal and the coion concentration will be lower inside the 
crystal. The higher counterion concentration and lower 
coion concentration will be such that the net charge den- 
sity in the crystal vanishes. If one can ignore specific 
ion interactions and the effects of non-ideality, the ion 
concentrations can be determined using Donnan's 1911 
theory, treating the protein-solution interface as a mem- 
brane. This was done in Ref. M, 

In the vicinity of the interface, the ion densities must 
vary smoothly between the interior of the crystal and the 
exterior salt solution. In particular, as shown in Fig. [I], 



one would expect counterions to spill out into the salt so- 
lution and coions to similarly spill into the protein, lead- 
ing to the formation of an electrical double layer across 
the interface. The space charge gives rise to an elec- 
tric field and a jump in the mean electrostatic potential 
between the salt solution and the protein crystal. This 
potential difference is sometimes termed the Donnan po- 
tential, or the Galvani potential jl2| . In Donnan's theory, 
this potential jump is self-consistently responsible for the 
change in ion concentrations in the bulk phases. 

In the theory of metals, a closely analogous phe- 
nomenon occurs at a metal surface where the electron 
density spills out into vacuum. The resulting double 
layer often makes an important contribution to the work 
function of metals, and was one of the early targets for 
density functional theory |l3| . In the metal case there is 
only one density, namely the electron density, which van- 
ishes in the vacuum. The electron density in the interior 
is very high and the electrons form a highly degenerate 
Fermi liquid, so that careful attention has to be paid 
to exchange and correlation effects. In contrast, in the 
protein case one has two densities, for counterions and 
coions, which approach a common value in the salt solu- 
tion outside the crystal. Moreover, for a 1:1 electrolyte in 
a high dielectric solvent, correlation effects are not that 
important and for simplicity in our calculation we ignore 
them. 

In the next section we determine the potential profile 
at the crystal-solution interface, and calculate the associ- 
ated free energy which represents the contribution of the 
electrical structure to the interfacial tension. The third 
section contains results of these calculations along with 
approximate analytic expressions for the magnitude of 
this contribution. The last section is a conclusion. 



II. CALCULATION 

We consider the case where we have a dense protein 
crystal coexisting with a very dilute protein solution, suf- 
ficiently dilute that the few protein molecules in it can be 
neglected and it treated as simply a salt solution. Rather 
than considering a specific protein crystal, we take over 
the jellium concept from the theory of metals to make a 
general estimate of the effect of salt on the interfacial ten- 
sion. We replace the detailed charge density due to the 
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FIG. 1: Ion density profiles for jellium model of a protein 
crystal-solution interface. The density profiles are computed 
from the numerical solution to the Poisson-Boltzmann equa- 
tion, and have been normalised by the salt concentration p a . 
In this example, the protein charge density p p (assumed pos- 
itive) was chosen equal to p s . The densities of counter and 
coions are p_ and p+, respectively. Distance is normalised by 
the Debye screening length defined in Eq. (H). 



protein molecules by a uniform background charge den- 
sity, cut off abruptly at the interface. Thus our model of 
the bare protein crystal comprises a uniformly charged 
half-space, with a charge density 
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Pp (* < 0) 
(z > 0) 



(1) 



where p p is the mean charge density in the protein crys- 
tal. For definitcncss we take the protein to be positively 
charged so the counterions are negative and the coions 
are positive; both are monovalent. The discontinuity in 
background charge density at z = is dressed by the 
counterion and coion densities P-(z) and p + (z). These 
satisfy p± — ► p a as Z — ► oo where p s is the salt concen- 
tration, and (p_ — p+) — > p p as z — > — oo. 

For notational convenience, we work in units where 
e = k^T = 1. In these units, the Coulomb potential 
energy U between a pair of elementary charges separated 
by r is U = Ib/t, where Ib is the Bjerrum length, equal to 
0.72 nm in water at room temperature (Ib = e /^neksT). 
We assume a constant value of Ib, and ignore dielectric 
effects. 

Below, we will find that in the limit of very large salt 
concentrations the potential is uniform everywhere, in 
the protein crystal and in solution, and so can be set 
equal to zero. In this limit the only relevant contributions 
to the free energy are independent of salt concentration. 
For the purposes of considering the effect of reducing the 
salt concentration, we assume that when the concentra- 
tion is reduced the protein crystal remains unchanged 
and interacts with the salt ions solely through a mean 
electrostatic potential, which we will calculate. We can 



then approximate the protein crystal by a uniform back- 
ground charge density, Eq. (0), and ignore the rest of the 
protein interactions until we come to estimate their con- 
tribution to the interfacial tension. These assumptions 
are reasonable if the interactions of the protein molecules 
are strong enough to form a dense, relatively rigid, crys- 
tal. 

We will use a grand potential, f2 e ; , which contains only 
ideal solution terms for the ions and the associated elec- 
tric field at the interface. We use a grand potential be- 
cause as usual the calculation of the interfacial tension 
is easiest at fixed chemical potential not density. The 
grand potential fl e i, is, 
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The first terms in u> are the ideal solution terms (the ions 
share a common chemical potential p± = log/? s ). The 
last term is the electrostatic energy, wherein E = —d<fi/dz 
is the electric field strength corresponding to an electro- 
static potential <f> which satisfies the Poisson equation, 
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(3) 



A correlation term Q co - [T [p±(z)] could also be included in 
Eq. (||) but, as already mentioned, for the purposes of the 
present calculation we omit correlation effects and work 
within the simple mean field theory. 

The variational principle SQ e i/S[p±(z)] = applied 
in this problem yields p±(z) = p s exp[=p0]. The elec- 
trostatic potential then satisfies the Poisson-Boltzmann 
equation, 



d 2 ^ 
dz 2 
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Equation (Q) can be integrated once with respect to <fi to 
obtain 



2k 2 (cosh cj) — cosh < 



dz 



-8-kIbPp (4> - 4>t>) (z<0) 



(6) 



2k 2 (cosh0 - 1) 



(z>0) 



where 4>d is the Donnan potential, the potential <fi in the 
bulk of the protein crystal, so <j) — > </>d as z — > — oo. 
The potential in the bulk of the solution is 0, <f> — > as 
z — > oo. 

For z > Eq. (|) can be integrated analytically Q 
to get <f> = 21og[(l + Ce- K ' z )/(l - Ce"^ 2 )] where C = 
tanh(0o/4) and </>o = (f>(0) is the value of the potential 
at the interface itself. For z < 0, it is not possible to 
integrate analytically, it must be done numerically. Now, 
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FIG. 2: Numerical factor in the interfacial tension, Eq. (|l7[), 
as determined from the full numerical solution to the Poisson- 
Boltzmann equation. 



a sharp dividing line between them. We therefore have 
to calculate the grand potential then subtract the grand 
potentials for the bulk protein and solution states, thus 

/oo 
dz [lj(z)-0(-z)u(-oo)-6(z)uj(+oo)] (12) 
-oo 

where 9{z) is the Heaviside step-function, and w(±oo) 
are the limiting values of the grand potential density at 
z — > ±oo. They may be obtained from the definition of 
w in Eq. @ and Eq. (@) for the crystal, 



uj(— oo) 
oj(+oo) 



-2p s (cosh fa 
-2p s 



i D sinh 4>v) 



(13) 



noting that in the bulk phases the electric field is zero. 
Now, using Eq. (|2|) together with p± = p s exp(=F0) and 
E(z) = —dfa'dz from Eqs. @ and (|io|), we can write co(z) 
as a function of <j). Using this together with Eq. ( fT3"| ) we 
can then write Eq. (12) as 



the charge densities for the co- and counterions in the 
bulk of the protein crystal (i.e., z — > — oo) are 



p±(-oo) = p s exp(=f<h)- 



(7) 



A 7el = 2p s 



dz fasmh (j> — 6(—z) sinh fa 



III. RESULTS 



(14) 



Charge neutrality imposes 



p+(-oo) + p P = p-(-oo). 



(8) 



Combining these 2 equations yields an equation for the 
Donnan potential 



sinh^D = rr~- 
2p s 



(9) 



and using this equation and Eq. (||) in Eq. (j^j we have 

(dfa'dz) 2 = 2k;;? [cosh - cosh</> D , . , . 

-(^-0 D )sinh0 D ] ^<U). (10) 



Matching dcf>/dz across z = using Eq. 
Eq. (O) for z < 0, gives 



for z > and 



<^o = + (1 — cosh D ) /sinh </> D . 



(11) 



This not only fixes the complete solution for z > 0, but 
can also be used as the starting point for a numerical 
integration of Eq. ( pTj| ) into the z < half-space. Note 
that if Eq. ( ^Tj| ) is used to determine d<p/dz, the negative 
root should be taken; also the term [. . .] in this equation 
has a geometric interpretation as the distance between 
cosh0 and its tangent at <f> = fa), and thus is always 
positive. 

The benefit of starting with an expression for the grand 
potential is that the electrical contribution to the inter- 
facial tension, A7 e l) is easily calculated. The interfacial 
tension is the difference between the actual grand poten- 
tial per unit area of the interface and that it would have 
if each of the 2 phases continued unperturbed right up to 



Before turning to the results of our full calculations, 
we discuss the high salt limit. If p s 3> p p then <j> <C 1 
everywhere and the Poisson-Boltzmann equation (Eq. (|^) 
for z > and Eq. @ for z < 0) can be linearised. The 
Donnan potential is fa — p p /2p s , the potential at z = 0, 
Eq. (|ll|), becomes fa — fa/2, and the solution to the 
linearised Poisson-Boltzmann equation is 



(<fc/2) [2 - exp( Ks z)] (z < 0) 
(fa)/2)exp(~K s z) (z>0) 



In Eq. (14) these yield for the interfacial tension 



A 7oI = -^ (Ps» Pp ). 



(15) 



(16) 



This already shows the main features of the full so- 
lution, being negative and of decreasing magnitude as 
the salt concentration increases. The fact that the result 
should be negative can be seen directly from Eq. jl^). 
If the ion density profiles were such that they remained 
constant right up to the interface where they had a sharp 
jump between the bulk values, then Eq. (12) would give 
= 0. The fact that the ion density profiles are re- 
laxed compared to such a trial density profile implies that 
A7oi < 0. Note that the electrostatic energy is always 
positive though (it is the last term in Eq. (0)). 

For the general case, we note that k~ is the only 
length scale that enters into the Poisson-Boltzmann equa- 
tion, thus dimensional analysis allows us to write 



2k s Pp 



,Ps 



(17) 
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FIG. 3: Predicted interfacial tension for lysozyme crystals, 
based on the present work and previous models ^, The 
dashed line is the high-salt limit 70 estimated in the text. The 
solid line incorporates the correction in Fig. ^. The dotted line 
is the result if / = 1 is used in Eq. (]17|). 



where f(x) is a numerical factor which amends the high 
salt scaling limit and depends only on the ratio x — 
p s / pp. We have implemented a numerical scheme to inte- 
grate Eqs. ([ll]) and (|l4|), and determine fix). The result 
is shown in Fig. |^. The function obviously obeys / — > 1 
as Ps/ P P — > 00, but the interesting thing is that / does 
not deviate greatly from unity even at much lower salt 
concentrations, for example / ps 1.4 at p s j p p ~ 0.1. This 
indicates that the high-salt scaling limit Eq. ([l6|) is not 
a bad approximation even at much lower salt concentra- 
tions, provided that we use the exact expression for the 
Donnan potential given in Eq. (j^). 

We can further approximate the contribution to the 
surface potential by using the high salt expression for 
(j>D, which is obtained by linearising Eq. (^J), and is 
0d = Pp/(2ps)- With this approximation we obtain 
A7oi = —pp/(8p s K s ). At high salt the contribution of 
the ions to the thermodynamic potential of the protein 
crystal is again obtained by linearisation and is given by 
uj{— 00) = — 2p s + p 2 /(4p s ). The first term is simply 
cj(+oo), it is linear in p s and has no effect, the second 
term is proportional to the square of the protein density 
and so is a, positive, contribution to the second virial 
coefficient. This has been noted before, see Refs. [|[ [l5[ . 
So, the contribution to the interfacial tension is of order 
of the contribution to the free energy density times the 
Debye length re J 1 , which defines the thickness of the in- 
terface of course. This is what we might have expected 
simply on the basis of dimensional analysis, although di- 
mensional analysis would not of course have told us that 
A7oi is negative. 

Our proteins are charged colloids and at least at high 
salt, charged colloids are often viewed as repelling each 
other via a Derjaguin-Landau-Verwey-Overbeek (DLVO) 
pairwise repulsion which has the form Q 2 Ib exp(—K 8 r)/r 



for a pair of colloids each with charge Q and r apart 
|]l4j| . Treating such repulsions within a mean- field ap- 
proximation yields a contribution to the free energy den- 
sity which is exactly what we have found, p 2 /(4p s ) [fl6|| . 
This free energy density is just a free energy per protein 
molecule of order Qp p / (Ap s ). If we assume an interface 
of width re J 1 at the edge of a crystal, within which there 
are (Pp/Q)^ 1 protein molecules per unit area of the in- 
terface, each of which has a DLVO energy which is some 
fraction of its value in the bulk, then the contribution of 
the DLVO repulsion to the interfacial tension is of order 
—p 2 /(p s n s ). Thus in the high salt limit our results for the 
contributions to the bulk free energy density and to the 
interfacial tension are of the same order as those obtained 
by assuming the protein molecules repel each other via 
a pairwise repulsion of the DLVO form; provided we use 
a mean-field approximation for the bulk and assume a 
protein free energy density profile of width re J 1 at the 
interface. Note that if we used a DLVO potential and 
assumed its contribution to A7 c i was simply due to a 
missing nearest-neighbour repulsion per molecule at the 
surface, then A7 c i = Q 2 Ib exp(— K s b)/b 3 , where b is the 
lattice spacing. This is rather different from our Donnan- 
potential result. Our finding that, in the high salt limit, 
our results obtained via a Donnan potential approach 
plus the Poisson-Boltzmann equation for the interface 
profile, can be obtained (up to numerical prefactors) via 
a DLVO potential is consistent with work on charged syn- 
thetic colloids. This work found that a DLVO pair poten- 
tial is adequate to describe the system, see for example 
Ref. By contrast, at low salt, charged colloids are 

not well described by a pair potential, see Refs. H II 
and references therein. 



A. Predictions for lysozyme 

To make these results more concrete, we now consider 
a specific example. We choose the lysozyme/NaCl/water 
system, since considerable information has been collected 
here. Elsewhere it has been shown that the crystallisa- 
tion boundary in this system and the second virial coef- 
ficient data can be fitted reasonably well with a model 
of hard spheres with sticky patches ||, |2^]. The sticky 
patches are short-ranged directional potential wells. Fol- 
lowing Ref. H we perform calculations for a well depth 
e = 7.4 fceT per patch. With attractions of this strength, 
and at high salt concentration, the protein crystal coex- 
ists with a very dilute protein solution. We estimate 
the interfacial tension in this model, in the high salt 
limit, to be 70 ps e/l 2 where I 2 is the interfacial area 
per protein pJ. Since the dimensions of lysozyme are 
4.5 x 3.0 x 3.0 nm 3 ps Z 3 , we estimate I 2 ps 12 nm 2 
and consequently 70 ps O.G/cB^nm" 2 = 2.5 mNm' 1 (by 
way of comparison, this is intermediate between the in- 
terfacial tension of a clean oil/water interface typically 
50mNm _1 , and interfacial tensions that have been mea- 
sured between colloidal phases in the range 1-20 /iNm -1 
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pl| . |22|]). To use our theory for the electrical correction to 
7o, we need the charge density in the crystal. Lysozyme 
has quite a high, pH-dependent charge, Q, and for the 
purpose of the present calculation we use Q = 10 |23f| 
which gives p p w Q/^ 3 « 0.25 nm -3 = 0.4 M. 

Fig. g shows the predicted surface tension, 7 = 70 + 
A7 c i, as a function of salt concentration. As might be 
expected, the electrical structure at the protein crystal- 
solution interface starts to significantly reduce the inter- 
facial tension for p s < 0.4 M. As the physiological salt 
concentration is 0.15 M, the interfacial tension of a pro- 
tein crystal in vivo is significantly reduced by the dou- 
ble layer. Our assumption that the density of protein 
molecules in the solution phase coexisting with the crys- 
tal is so low that it can be approximated by zero, is rea- 
sonable down to a salt concentration of about 0.1 M, with 
sticky patches of strength e « 7.4 k^,T ||. Below this we 
would need to take account of the density of protein in 
the solution phase coexisting with the crystal. 

We have been unable to find any experimental determi- 
nation of the protein crystal-solution interfacial tension, 
but we hope that our calculations may stimulate exper- 
imental work to confirm our results. However, a predic- 
tion can be extracted from our theory for the variation 
in protein solubility, i.e., the lysozyme concentration in 
the fluid phase coexisting with the crystal, as a function 
of the concentration of NaCl (if], [24J. As noted above, 
at high salt concentration the grand potential density in 
the protein crystal is u>{— 00) = — 2p s + pi/(Ap s ). This 
implies a contribution to the excess chemical potential of 
a protein molecule of Qp p /(2p s ). Treating the solution 
phase which coexists with the crystal as ideal (since the 
protein concentration is very small) this variation of the 
excess chemical potential with salt concentration yields 
for the solubility as a function of salt concentration 



lnp sol (p s ) = lnp so i(p s -> 00) + 



Qp P 

2p s 



(18) 



where p so i(ps) is the charge density of the protein 
molecules in the dilute solution which coexists at equi- 
librium with the protein crystal. This recovers a result 
also found for the high salt limit of the theory in Ref . || . 
As the first term on the right-hand side is a constant, 
our prediction is that the logarithm of the solubility is 
inversely proportional to the salt concentration. This 
is same dependence as found in experiment [24| but, as 
noted in Ref. ||, the slope predicted by Eq. (18) is close 
to twice the slope measured for the experimental data. 



B. Comparison with surface-force apparatus 
measurements 

A closely related problem to that of the free energy of a 
single crystal surface is that of the potential of mean force 
per unit area W(s) between a pair of protein surfaces as 
their separation s is varied. This is related to the inter- 
facial tension problem since W(0) = —2j if we choose 



the reference state such that W(oo) = 0. Its derivative 
with respect to separation is the negative of the force per 
unit area between two surfaces. We expect that the im- 
portant contribution to W from the sticky patches will 
be operative only when s < 5, where S is the range of 
the sticky patch potential. Previous work on the phase 
diagram suggests S is only a few percent of the protein 
diameter || [2CJ]. If s > S, W arises from the overlap of 
the double layers and should be repulsive since both sur- 
faces carry the same charge. In fact, we would expect W 
to exhibit a potential barrier at s ~ 5, of height — 2A7 i 
and extending out to a distance of order kJ 1 . 

The interaction between protein surfaces is more ex- 
perimentally accessible than the interfacial tension, as it 
can measured via a direct force experiment. Sivasankar 
et al have recently looked at the interaction between 
streptavidin-covered surfaces using a surface force appa- 
ratus fzsj] . This is not quite our system as we consider a 
bulk crystal whereas the experiments are on streptavidin 
monolayers, but so long as the Debye screening length 
is smaller than the protein diameter, the thickness of 
the layer over which the ion densities vary will not be 
much larger than the thickness of a monolayer of protein 
molecules. Sivasankar et al found that the interaction 
was well described by the non-linear Poisson-Boltzmann 
equation with a fixed effective surface charge density of 
order 10mCm~ 2 . Whilst the present theory does not 
calculate the interaction between a pair of surfaces, we 
note that the electrostatic potential outside the crystal is 
the same as would be obtained for a surface charge den- 
sity a — (4p s / k s ) sinh(0o/2), where <po is the potential at 
the surface determined earlier (the simplest way to derive 
this result is via matching the electric field strength at 
z = 0). For example, for lysozyme at p s = p p = 0.4 M, 
we find, using this equation, an effective surface charge 
density a « 0.06enm~ 2 w 10mCm~ 2 , of exactly the 
same magnitude seen by Sivasankar et al j2S| . 

Finally we comment briefly on the case of a single pro- 
tein molecule approaching the surface of a protein crystal. 
In this case one would expect a barrier of height Q4>q 
to be present in the potential of mean force, before the 
isolated protein encounters the short range attraction. 



IV. CONCLUSION 

We have estimated the contribution to the interfacial 
tension of a protein crystal made by the electrical double 
layer at the interface. We have used a jellium model in 
which the protein is replaced by a uniform background 
charge density, and have solved the Poisson-Boltzmann 
equation for the ion density profiles for a jellium half- 
space. Whilst a full solution cannot be obtained analyt- 
ically, we find a good approximation is to use the high- 
salt limit result in Eq. (|l^) with the Donnan potential 
given in Eq. (^|). The contribution of the double layer 
is, perhaps surprisingly, negative and for a protein not 
too close to its isoelectric point and of a typical size, it 



6 



amounts to roughly 10% of the total interfacial tension at 
salt concentrations of order 0.1 M, see Fig. ||. As the salt 
concentration is increased the contribution of the double 
layer decreases, due to the decreasing step in potential 
as the interface is crossed and its decreasing width. We 
have also briefly considered the potential of mean force 
between a pair of protein surfaces as they are brought to- 
gether. The mean force is repulsive and agrees with that 
found by Sivasankar et al pq ] for the repulsion between 
monolayers of the protein streptavidin. 

The fact that the contribution of the double layer 
to the interfacial tension is negative and largest at low 
salt concentrations suggests that proteins will crystallise 
more readily at low rather than at high concentrations 
of a salt such as NaCl. At lower salt concentrations the 
interfacial tension will be lower, all other things being 
equal, and the rate of nucleation varies as the exponen- 
tial of minus the cube of the interfacial tension. However, 
decreasing the concentration of salt will make the protein 
crystal more soluble, it stabilises the solution phase at 
the expense of the crystal || . Thus if the salt concentra- 
tion is too low the solution will be the equilibrium phase 
and crystallisation will not then occur at all. Because 
lowering the salt concentration lowers both the relative 
stability of the crystal with respect to the solution and 
the interfacial tension, we conclude that there is no clear 
'best' salt concentration at which to attempt to crys- 



tallise a protein. Also, note that by salt we mean a salt 
such as NaCl in the lysozyme system, in which there are 
no specific ion effects such as binding of an ion to the 
protein. 

As well as the interaction between a pair of protein 
surfaces, we mention some other directions in which our 
theory could be extended. Firstly, ion correlations and 
the effects of an inhomogeneous charge distribution could 
be included. Secondly, we have neglected excluded vol- 
ume effects and changes in dielectric properties which 
are potentially important in the interior of the protein 
crystal. Thirdly, we have only considered the case where 
the protein crystal coexists with a dilute protein solu- 
tion, such that the protein concentration in the solution 
makes a negligible contribution to the charge balance. It 
has been suggested that protein crystallisation is facili- 
tated by a metastable liquid-liquid demixing transition 
jl], [H , in which case one should certainly examine the 
effect of a more concentrated protein solution phase ||. 

We would like to dedicate this paper to Jean-Pierre 
Hansen, on the occasion of his 60th birthday. Indeed, 
this is particularly appropriate since it is a pleasure to 
record that the present calculations were started during 
the recent conference to celebrate Jean-Pierre's many sci- 
entific contributions to liquid state theory. The work of 
RPS is supported by EPSRC (GR/N36981). 
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